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Abstract. A model for traffic flow in street networks or material flows in 
supply networks is presented, that takes into account the conservation of cars 
or materials and other significant features of traffic flows such as jam formation, 
spillovers, and load-dependent transportation times. Furthermore, conflicts or 
coordination problems of intersecting or merging flows are considered as well. 
Making assumptions regarding the permeability of the intersection as a function 
of the conflicting flows and the queue lengths, we find self-organized oscillations 
in the flows similar to the operation of traffic lights. 

1. Introduction. Material flows are found in many places of the world. This con- 
cerns, for example, traffic flows in urban areas or flows of commodities in logistic 
systems. There is also some similarity with material flows in production or bio- 
logical systems, from cells over bodies upto ecological food chains. Many of these 
material flows are not of diffusive nature or going on in continuous space. They 
are often directed and organized in networks. In comparison with data flows in 
information networks, however, there are conservation laws, which can be used to 
set up equations for material flows in networks. It turns out, however, that this 
is not a trivial task. While there is already a controversial discussion about the 
correct equations representing traffic flows along road sections [51 [TTJ [33J j their 
combination in often complex and irregular networks poses further challenges. In 
particular, there have been several publications on the treatment of the bound- 
ary conditions at nodes (connections) of several network links (i.e. road sections) 
[251[3[2i[Tl[T0l[S[l[21[23j[24i- In particular, the modelling of merging and 
intersecting flows is not unique, as there are many possible forms of organization, 
including the use of traffic lights. Then, however, the question comes up how these 
traffic lights should be operated, coordinated, and optimized. In order to address 
these questions, in Sec. [2] we formulate a simple model for network flows, which 
contains the main ingredients of material or traffic flows. Section [3J will then discuss 
the treatment of diverges, merges, and intersections. Equations for the interaction- 
dependent permeability at merging zones and intersections will be formulated in 
Sec. HJ We will see that, under certain conditions, they lead to spontaneous oscilla- 
tions, which have features similar to the operation of traffic lights. Finally, Sec. [5] 
summarizes and concludes this paper. 
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2. Flows in Networks. The following section will start with a summary of the 
equations derived for traffic flows in networks in a previous paper. These equations 
are based on the following assumptions: 

• The road network can be decomposed into road sections of homogeneous ca- 
pacity (links) and nodes describing their connections. 

• The traffic dynamics along the links is sufficiently well described by the 
Lighthill-Whitham model, i.e. the continuity equation for vehicle conserva- 
tion and a flow-density relationship ("fundamental diagram"). This assumes 
adiabatic speed adjustments, i.e. that acceleration and deceleration times can 
be neglected. 

• The parameters of vehicles such as the maximum speed V® and the safe time 
headway T are assumed to be identical in the same road section, and who 
enters a road section first exits first (FIFO principle). That is, overtaking is 
assumed to be negligible. 

• The fundamental diagram can be well approximated by a triangular shape, 
with an increasing slope V® at low densities and a decreasing slope c in the 
congested regime. This implies two constant characteristic speeds: While V® 
corresponds to the free speed or speed limit on road section i, 



the dissolution speed of the downstream front of a traffic jam and the velocity 
of upstream propagation of perturbations in congested traffic. While p max 
denotes the maximum vehicle density in vehicle queues, T ~ 1.8s is the safe 
time gap between two successive vehicles. 
• The vehicle density in traffic jams is basically constant. 
These assumptions may be compensated for by suitable corrections [llj , but already 
the model below displays a rich spectrum of spatio-temporal behaviors and contains 
the main elements of traffic dynamics we are interested in here. 

2.1. Flow Conservation Laws. In the following, we will introduce our equations 
for traffic flows in networks very shortly, as a detailed justification and derivation 
has been given elsewhere [12j [T3l HI]- These equations are also meaningful for 
pipeline networks [3] (if complemented by equations for momentum conservation), 
logistic systems [21], or supply networks [19] . Our notation is illustrated in Fig. [I] 
Compared to Ref. [14] . we will use a simplified notation, here0 The arrival flow 
Aj(t) denotes the actual inflow of vehicles into the upstream end of road section j, 
while Oj(t) is the actual departure flow, i.e. the flow of vehicles leaving road section 
j at its downstream end. The quantity 



represents the maximum in- or outflow of road section j. All the above quantities 
refer to flows per lane. Ij is the number of lanes and Lj the length of road section 
j. lj(t) < Lj is the length of the congested area on link j (measured from the 
downstream end), and ANj is the number of stopped or delayed vehicles, see Eqs. 

The arrival flow Aj (t ) has previously been denoted by Q| rr (t ) , the potential arrival flow A j (t ) 
by Q^ 11,pot (t), the departure flow Oj(t) by Q^ cp (t) and the potential departure flow Oj(t) by 
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(fT6|) and (fl4|) . With these definitions, we can formulate constraints for the actual 
arrival and departure flows, which are given by the potential arrival flows Aj(t) and 
the potential departure flows Oi(t), respectively. 

The actual arrival flow Aj (t) is limited by the maximum inflow Qj , if road section 
j is not fully congested (lj(t) < Lj). Otherwise (if lj = Lj) it is limited by the actual 
departure flow Oj(t — Lj/c) a time period Lj/c before, as it requires this time period 
until the downstream flow value has propagated upto the upstream end of the road 
section by forward movement of vehicles under congested traffic conditions. This 
implies 

0<A jW <l jW :={| (< _ Vc) (3) 

Moreover, the potential departure flow Oi(t) of road section i is given by its perme- 
ability 7i(t) times the maximum outflow Qi from this road section, if vehicles are 
queued up (AiVj > 0) and waiting to leave. Otherwise (if AiV$ = 0) the outflow is 
limited by the permeability times the arrival flow Ai a time period Li/V® before, as 
this is the time period that entering vehicles need to reach the end of road section 
i when moving freely at the speed V®. This gives the additional relationship 

< Oi(t) < Oi(t) .- 7 i(t) | Q, if AN , {t) >Q (4) 

The permeability ji(t) for traffic flows at the downstream end of section i can 
assume values between and 1. In case of a traffic light, ~fi(t) = 1 corresponds to a 
green light for road section i, while "fi(t) = corresponds to a red or amber light. 
Alternatively and shorter than Eqs. {3]) and ((4|) one can write 

A, (t) = max O (lj (t) < Lj), Oj (t - Lj /c)] (5) 

and 

dt(t) = fi(t) max [QSiANi > 0), A t {t - L,/^ )] , (6) 

where the Heaviside function O is 1, if the argument (inequality) has the logical 
value "true" , otherwise it is 0. Note that the above treatment of the traffic flow in a 
road section requires the specification of the boundary conditions only, as we have 
integrated up Lighthill's and Whitham's partial differential equation over the length 
of the road section. The dynamics in the inner part of the section can be easily 
reconstructed from the boundary conditions thanks to the constant characteristic 
speeds. However, a certain point of the road section may be determined either 
from the upstream boundary (in the case of free traffic) or by the downstream 
boundary (if lying in the congested area, i.e. behind the upstream congestion front). 
Therefore, we have a switching between the influence of the upstream and the 
downstream boundary conditions, which makes the dynamics both, complicated 
and interesting. This switching results from the maximum functions above and 
implies also that material flows in networks are described by hybrid equations. 
Although the dynamics is determined by linear ordinary differential equations in 
all regimes, the switching between the regimes can imply a complex dynamics and 
even deterministic chaos [32) . 

Complementary to the above equations, we have now to specify the constraints 
for the nodes, i.e. the connection, merging, diverging or intersection points of the 
homogeneous road sections. Let the ingoing links be denoted by the index i and the 
outgoing ones by j. To distinguish quantities more easily when we insert concrete 
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values 1,2,... for i and j, we mark quantities of outgoing links additionally by a 
prime ('). 

Due to the condition of flow conservation, the arrival flow into a road section j 
with Ij lanes must agree with the sum of the fractions onj of all outflows ijOj(i) 
turning into road section i. Additionally, the arrival flows are limited, i.e. we have 

r A (*) = J2 w*Hj ^ l 'A (*) ( 7 ) 

i 

for all j. Of course, the turning fractions ay > are normalized due to flow 
conservation: 

EM*) = 1 - (8) 

j 

In cases of no merging flows, Eq. ([7]) simplifies to 

Jj^ (t) - IiOi{t) aij < I'jA'j (t) (9) 

for all j. At the same time, < Oj(t) < Oi(i) must be fulfilled for all i. Together, 
this implies 



Oi(t) < min 



Oi{t), min 



(10) 



for all z. 

The advantage of the above model is that it contains the most important ele- 
ments of the traffic dynamics in networks. This includes the transition from free to 
congested traffic flows due to lack of capacity, the propagation speeds of vehicles and 
congested traffic, spillover effects (i.e. obstructions when entering fully congested 
road sections) and, implicitly, load-dependent travel times as well. 

2.2. Two Views on Traffic Jams. Let us study the traffic dynamics on the road 
sections in more detail. Traffic jams can be handled in two different ways: First by 
determining the number of cars that are delayed compared to free traffic or, second, 
by determining fronts and ends of traffic jams. The former method is more simple, 
but it cannot deal correctly with spill-over effects, when the end of a traffic jam 
reaches the end of a road section. Therefore, the first method is sufficient only in 
situations where the spatial capacity of road sections is never exceeded. 

2.2.1. Method 1: Number of Delayed Vehicles. The first method just determines 
the difference between the number N\ n (t) of vehicles that would reach the end of 
road section i upto time t and the number N° ut (t) of vehicles that actually leave 
the road section upto this time. Nl n (t) just corresponds to the number of vehicles 
which have entered the road section upto time t — Li/V®, as Li/V® is the free travel 
time. This implies 



Nt(t) = J dt> Mt> -Li/V?), (11) 
o 

while the number of vehicles that have acually left the road section upto time t is 

t 

N° ut (t)= I dt 1 Oi(t'). (12) 
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Hence, the number ANi(t) of delayed vehicles is given by 

t 

ANi(t) = J dt' [Mt' - Li/V?) - Oi(t')} > . (13) 
o 

Alternatively, one can use the following differential equation for the temporal change 
in the number of delayed vehicles: 

dAN< 



Ai(t-Li/V?)-Oi(t). (14) 

In contrast, the number of all vehicles on road section i (independently of whether 
they are delayed or not) changes in time according to 

< ^ = A i {t)-O i (t). (15) 

2.2.2. Method 2: Jam Formation and Resolution. In our simple macroscopic traffic 
model, the formation and resolution of traffic jams is described by the shock wave 
equations, where we have the two characteristic speeds V® (the free speed) and c 
(the jam resolution speed). According to the theory of shock waves [291 |3"T] . the 
upstream end of a traffic jam, which is located at a place U{t) > upstream of the 
end of road section i, is moving at the speed 

dk = Mt- [Li - k(t)]/V-) -Oi{t- W)/c) 
dt Pl (t)-p 2 (t) [ ' 

with the (free) density 

p 1 (t)=A i (t-[L i -l i (t)]/V°)/V° (17) 

immediately before the upstream shock front and the (congested) density 

P2 (t) = [1 - TOi(t - k(t)/c)] Piaax (18) 

immediately downstream of it. This is, because free traffic is upstream of the shock 
front, and congested traffic downstream of it (for details see Eqs. (1.6) and (1.4) in 
Ref. [33]). In contrast, the downstream front of a traffic jam is moving at the speed 

0-Oi(t-U(t)/c) Oi(t-k(t)/c) 



Oi(t - k(t)/c)/VP Pmax -Oi(t- li(t)/c)/V? ' 



(19) 



since congested traffic with zero flow is upstream of the shock front and free traffic 
flow occurs downstream of it. 

2.2.3. Comparison of the Two Methods. Let us discuss a simple example to make 
the differences of both descriptions clearer. For this, we assume that, at time t = 0, 
traffic flow on the overall road section i is free, i.e. any traffic jam has resolved and 
there are no delayed vehicles. The flow shall be stopped by a red traffic light for 
a time period t$. At time t = to, the traffic light shall turn green, and the formed 
traffic jam shall resolve. For the arrival flow, we simply assume a constant value 
Ai, and the road section shall be long enough to take up the forming traffic jam. 
Moreover, the departure flow shall be Oi. Then, according to method 1, the number 
of delayed vehicles at time t is 

AJVi(*o) = Ait , (20) 

and it is reduced according to 

AN i (t) = A i t Q -(O i -A i )(t-t ). (21) 



6 



DIRK HELBING, JAN SIEGMEIER, AND STEFAN LAMMER 



Therefore, any delays are resolved after a time period 

t - t Q = Att ° = ANl{h) , (22) 

i.e. at time 

h = t0 O^Al- (23) 

Afterwards, AN,(t) = 0. 

In contrast, the end of the traffic jam grows with the speed 

dk = A-o = l = 

dt ~ Ai/V? - (1 - 0) Pmax Pmax/A, - 1/V° ' 1 " [ ' 

Therefore, we have k(to) = C^o- Surprisingly, this is greater than ANi(to) / p max , 
i.e. the expected length of the traffic jam based on the number of delayed vehicles. 
The reason is that the delay of a vehicle joining the traffic jam at location Xi = Li — k 
is noticed at the downstream end of the road section only after a time period h/V®. 
The resolution of the traffic jam starts from the downstream end with the speed 

0-Qi -1 



(25) 



Pmax - Qi/VP Pmax /Qi - l/V? 

if the outflow is free (i.e. Oi = Qi), otherwise with the speed 

— 7 _ n I \ = ~ c ' ( 26 > 

Pmax \Pmax ^i/C) 

since congested traffic with zero flow and maximum density is upstream of the shock 
front. 

Obviously, the jam resolution has reached the further growing, upstream jam 
front when Cit = c(t — to). Therefore, the jam of density p max has disappeared 
after a time period t — to = Cito/{c — C'i), i.e. at time 

h = ct /(c - Ci) . (27) 

Surprisingly, it can be shown that t\ < t 2 , i.e. the traffic jam resolves before the 
number of delayed vehicles reaches a value of zero. In fact, it still takes the time 
Crfi/V® until the last delayed vehicle has left the road section, where 

MIV? - (Pmax - OJc) 

is the shock front between free upstream traffic flow and the congested outflow Oi, 
which usually differs from the speed Vi = Oj/[(l — TQi)p max ] of outflowing vehicles. 
For O t = Q t , we have V® = V? because of 1/c = p max /Qi - l/V®. 

Undelayed traffic starts when this shock front reaches the end of the road section, 
i.e. at time 

„ -Jl + °L\- r^rr ( 1 + Ct(A ' /V - ^ fe=i + Ct0 ' /C ) ■ (29) 

^ VP J 1 - C l /c V Ai - O t ) 

Inserting d(Ai/VP — p max ) — —Ai eventually gives t 2 = t Oi/Oi — Ai). This agrees 
perfectly with the above result for the first method (based on vehicle delays rather 
than traffic jams). 

In conclusion, both methods of dealing with traffic jams are consistent, and 
delayed vehicles occur as soon as traffic jam formation begins. However, according 
to method 1, a queued vehicle at position Xi — Li — k is counted as delayed only 



^ = 7 7T70 „ 73 (28) 
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after an extra time period h/V® , but it is counted as undelayed after the same extra 
time period. This is because method 1 counts on the basis of vehicle arrivals at the 
downstream end of road section i. 

As it is much simpler to use the method 1 based on determining the number of 
delayed vehicles than using method 2 based on determining the movement of shock 
fronts, we will use method 1 in the following. More specifically, in Eq. (|3|) we will 
replace lj(t) < L, by ANj(t) < Af ax := L jPmax and l 3 (t) = L. } by ANj(t) = N™ a *. 
This corresponds to a situation in which the vehicles would not queue up along the 
road section, but at the downstream end of the road section, like in a wide parking 
lot or on top of each other. As long as road section j is not fully congested, 
this difference does not matter significantly. If it is fully congested, the dynamics 
will potentially be different, defining a modified model of material network flows. 
However both, the original and the modified model fulfill the conservation equation 
and show spillover effects. 



2.2.4. Calculation of Cumulative and Maximum Individual Waiting Times. In Ref. 
[T2] . we have derived a delay differential equation to determine the travel time Tj(i) 
of a vehicle entering road section i at time t (see also Ref. [IJ HI IS]): 

dTjjt) Mi) _ 

dt Oi(t + Ti(t)) ■ [ 1 

According to this, the travel time Ti(t) increases with time, when the arrival rate 
Ai at the time t of entry exceeds the departure rate Oi at the leaving time t + Ti(t), 
while it decreases when it is lower. It is remarkable that this formula does not 
explicitly depend on the velocities on the road section, but only on the arrival and 
departure rates. 

Another method to determine the travel times is to integrate up over the number 
of vehicles arriving in road section i, 

i 

N t A (t) = J dt' Mt') = Nt(t + U/V?) , (31) 
o 

and over the number of vehicles leaving it, 

t 

N°(t)= [dt'O i (t') = Nr t (t), (32) 



starting at at time t = when there are no vehicles in the road. If T!(t) denotes 
the time at which N?(t + Tftt)) = N^(t), then T/(t) is the travel time of a vehicle 
entering road section i at time t and 

T i (t)=T!(t)-L i /V? (33) 



is its waiting time. 
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Another interesting quantity is the cumulative waiting time T?(t), which is de- 
termined by integrating up over the number AiVj of all delayed vehicles. We obtain 



t 



Tf{t) = J dt' ANi(t') = J dt' [Nl n (t' - Li/Vt) - N° ut {t')} 

o o 
t t' 

dt' J dt" [Ai{t" - Li/Vi) - 00')} (34) 
o o 

and the differential equation 

t 

= ANi(t) = J dt' [Mt 1 - Li/V°) - O l (t')} . (35) 
o 

For a constant arrival flow A; and a red traffic light from t — to t — to (i.e. 
Oi{t) =0), we find 

At 2 

T^^f. (36) 

In this time period, a number of Ni(t) = Ait vehicles accumulates, which gives an 
average waiting time of 

AN(t a ) - 2 {67) 

at the end of the red light. The first vehicle has to wait twice as long, namely, a 
time period to. 



3. Treatment of Merging, Diverging and Intersection Points. While the 
last section has given general formulas that must be fulfilled at nodes connecting 
two or more different links, in the following we will give some concrete examples, 
how to deal with standard elements of street networks. For previous treatments of 
traffic flows at intersections see, for example, Refs. [51 [231 [71 [25], 



(b) 



(c) 
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Figure 1. Schematic illustration of the (a) diverging, (b) merging, 
and (c) intersecting flows discussed in this paper. 
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3.1. Diverging Flows: One Inflow and Several Outflows. In the case of one 

road section i diverging into several road sections j (see Fig. [T^,), Eqs. (flU)) and (JSJ) 
to imply 

U 



Oi (t) < min < 7^ (t) max 



Qi&(ANi > 0),Ai t 



V° 



mm 

j 



1% 0£%j 



max 



(QjQQjKL^Ojit-Lj/c) 



(38) 



for all i. If we assume that downstream road sections are never completely con- 
gested, this simplifies to 



Oi(t) = min |Qj, 



"fi max 



with 



Otherwise 



Qi = min 

3 



Q i e(AN i >0),A i (t-£i/t?)]} 



IiOLij 



mm 

j 



e(i 3 < Lj), 



I I QLij 



(39) 
(40) 

(41) 



3.2. Merging Flows: Two Inflows and One Outflow. We assume a flow 
l\0\(t) that splits into two flows I\Oi{t)an (going straight) and I\Oi(t)ai2 (turn- 
ing right), but a right-turning flow -^(^(i) merging with flow IiOi(t)an, as in 
turn-right-on-red setups (see Fig.IT})). For this situation, we have the equations 



I'Mt) = I l O l {t)a ll +I 2 2 (t)<I' 1 A' 1 {t), 
I' 2 A' 2 {t) = JiOi(t)ai 2 <I 2 A' 2 (t). 



One can derive 



and 



Let us set 



< Oi = min 



0<O 2 = 



Oi (*), 



I{A[(t) ~ l 2 2 (t) I' 2 A' 2 (t) 



2 (t) 



I%an hai 2 
I'Mt) - hOi{i)au- 



Oi 



Oi 



I[A[(t) I' 2 A' 2 {t) 



and 



2 {Oi) = min 



o 2 (t), 



I[A[(t) - hOiau 



(42) 
(43) 

(44) 

(45) 
(46) 
(47) 



Then, it can be shown that 2 (t) > and Oi(t) < [I{A[(t) - l 2 2 (t)}/{han), as 
demanded. If Oi (t) is chosen a value AOi smaller than specified in Eq. ([46]) , but 
2 is still set to the maximum related value 2 (Oi — AOi) according to Eq. (|47|) . 
the overall flow 

F = hOi+I 2 2 (48) 

is reduced as long as an < 1, since this goes along with additional turning flows 
(while the number of lanes does not matter!). Therefore, it is optimal to give 
priority to the outflow Oi(t) according to Eq. (|4*6]) and to add as much outflow 
2 (t) as capacity allows. This requires suitable flow control measures, otherwise 
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the optimum value of the overall flow F could not be reached. In fact, the merging 
flow would "steel" some of the capacity reserved for the "main" flow (i — 1), which 
would reduce the possible outflow Oi (t) and potentially cause a breakdown of free 
traffic flow, as it is known from on-ramp areas of freeways [34j . 

3.3. A Side Road Merging with a Main Road. Compared to the last section, 
the situation simplifies, if we have just a side road or secondary turning flow merging 
with a the flow of a main road without any turning flow away from the main road. 
In this case, we have an = 1 and a\ 2 = 0, which leaves us with the relationships 



0\ = min 

and 



Ox, 



i'Mt) 



h 



(49) 



0,(0,) = mm \d 2 (t), VMt \ h0l \ . (50) 

L 1 2 J 

according to Eqs. |gBJ and ((47)) . 

3.4. Intersection-Free Designs of Road Networks. With the formulas for the 
treatment of merges and diverges in the previous sections, it is already possible 
to simulate intersection-free designs of urban road networks, which do not need 
any traffic light control. The most well-known design of intersection-free nodes are 
roundabouts (see the upper left illustration in Fig. [2]). It is, however, also possible to 
construct other intersection-free designs based on subsequent merges and diverges 
of flows with different destinations. Two examples are presented in Fig. (2b and c. 

Although intersection-free designs require the driver to take small detours, such 
a road network will normally save travel time and fuel, given that the traffic volume 
is not too low. This is because intersections then need to be signalized in order to be 
safe and efficient H Traffic signals, however, imply that vehicles will often be stopped 
for considerable time intervals. This causes significant delays, at least for vehicles 
not being served by a green wave. Intersection-free designs, in contrast, do not 
necessarily require vehicles to stop. Therefore, the average speeds are expected to 
be higher and the travel times lower than for road networks with intersections. This 
has significant implications for urban transport planning, if intersections cannot be 
avoided by bridges or tunnels. 

3.5. Two Inflows and Two Outflows. The treatment of intersecting flows is 
more complicated than the treatment of merges and diverges. Moreover, the result- 
ing flows are only uniquely defined, if additional rules are introduced such as the 
optimization of the overall flow. Let us here treat the case of an intersection with 
two inflows and two outflows (see Fig. [TJ;). Equation implies the inequalities 

< I[A[(t) = hO^an + l 2 2 (t)a 21 < I'^t) , 

< I 2 A 2 (t) = h0 1 (t)a 12 +I 2 2 {t)a 22 < I 2 A' 2 {t) (51) 

with the constraints 

0<Or(t) < 6i(t), 

0< 2 {t) < 2 (t), (52) 



2 Of course, a first-come-first-serve or right-before-left rule will be sufficient at small traffic 
volumes. 
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Figure 2. Three examples for intersection- free designs of urban 
road networks. 



so that IjA'j(t) > is automatically fulfilled. The constraints (|52|) define an rectan- 
gular area of possible CVvalues in the 0\-Oi plane, where the size of the rectangle 
varies due to the time-dependence of Oi(t). The inequalities (fSTj) can be rewritten 
as 



and 



2 (t)< I ^- h °^ =:a 2 - b2 O l (t). 



(53) 



(54) 



They potentially cut away parts of this rectangle, and the remaining part defines 
the convex set of feasible points (0\, O2) at time t. We are interested to identify 
the "optimal" solution (01,0^), which maximizes the overall flow 



E 



/;4(t)=5>°*(*) 



(55) 



As this defines a linear optimization problem, the optimal solution corresponds to 
one of the corners of the convex set of feasible points, namely the one which is 
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touched first by the line 

Z - IjPi 

2 = 7 , (56) 

H 

when we reduce Z from high to low values. 

Let us, therefore, determine all possible corners of the convex set and the con- 
ditions, under which they correspond to the optimal solution. We will distinguish 
the following cases: 

(a) None of the boundary lines (|53|) and ([54)) corresponding to the equality signs 
cuts the rectangle defined by < O x (t) < 6 x (t) and < 2 (t) < 2 (t) in 
more than 1 point. This case applies, if a\ — b\0\ > 2 and a 2 — b 2 Oi > 2 , as 
a; > and hi > implies that both lines are falling or at least not increasing. 
Since the line (|56|) reflecting the goal function is falling as well, the optimal 
point is 

(0* 1 ,0* 2 ) = (0 1 ,0 2 ), (57) 

i.e. the outer corner of the rectangle corresponding to the potential or maxi- 
mum possible departure flows (see Fig. [3|. 




Figure 3. Illustration of the possible optimal solutions for two 
intersecting flows (see text for details). 



(b) Only one of the two boundary lines/border lines, 2 (t) — a\ — b\0\(t) or 
2 (t) = a 2 — b 2 Oi(t), cuts the rectangle in more than one point. Let us 
assume, this holds for line i, i.e. — biO\ < 2 . Then, the left cutting point 

(of o 1 ]) = { _ d 2 )/bi,d 2 ^j if ai > 6 2 , ^ 58 ^ 

I (0, dj) otherwise 

is the optimal point if Ii/I 2 < i.e. if the slope Ii/I 2 of the goal function 
(|56p is smaller than the one of the cutting border line. Otherwise, \ili/l 2 > 6j, 
the optimal point is given by the right cutting point 

(Of , 0£) = { S\ a * T bldl } if ? l >hdl ' (59) 
I (a»A,0) otherwise v ; 



(see Fig. [3]). 
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(c) If both border lines cut through the rectangle, but one of them lies above the 
other line, then only the lower line determines the optimal solution, which can 
be obtained as in case (b). Case (c) occurs if 0,2 — 62 O \ l > a\ — bi O J 1 and 
ei2 — b20\ r > ai — biO\ T (line 1 is the lower one) or if 02 — 62OJ 1 < at — biO{ 1 
and 02 — b20\ r < a\ — biO{ r (line 2 is the lower one). 

(d) The boundary lines cut each other in the inner part of the rectangle. This 
occurs if a\ — b\0\ < O2 and 02 — &2O1 < C>2- Then, the left-most cutting 
point (Of, 0\ ) is the optimal solution, if the slope I1/I2 of the goal function 
is smaller than the smallest slope of the two boundary lines, while it is the 
lower right cutting point (Of ,0 2 r ), if /1//2 is greater than the steepest slope 
of the two boundary lines, otherwise, the cutting point of the two boundary 
lines, 

' a 2 - ai aib 2 — bia 2 " 



(o[,o' 2i 



bi ' 62 — b\ 



(60) 



is the optimal point (see Fig. [3]). Mathematically speaking, we have 



(01,0*2) = I 



(o[,o> 2 ) 

(Of, Of) 
(Of, Of) 

(o[,o> 2 ) 
(o\\of) 



if h/I 2 < h < b 2 , 
if 61 < h/I 2 < b 2 , 
if 61 < b 2 < h/h, 
if h/h <b 2 < h, 
if 6 2 < h/h < bi, 
if 62 < bi < h/h, 



(61) 



It is astonishing that the simple problem of two intersecting traffic flows has 
so many different optimal solutions, which sensitively depend on the parameter 
values. This can reach from situations where both outgoing road sections experience 
the maximum possible outflows upto situations, where the outflow in the system- 
optimal point becomes zero for one of the road sections. A transition from one 
optimal solution to another one could easily be triggered by changes in the turning 
fractions entering the parameters at and bi, for example due to time-dependent 
turning fractions a.y(t). 



3.6. Inefficiencies due to Coordination Problems. An interesting question is 
how to actually establish the flows corresponding to the system optima that were 
determined in the previous sections on merging and intersecting flows. Of course, 
zero flows can be enforced by a red traffic light, while maximum possible flows can 
be established by a node design giving the right of way to one road (the "main" 
road). However, it is not so easy to support an optimimum point corresponding to 
mixed flows, such as (0^,0 2 )- That would need quite tricky intersection designs 
or the implementation of an intelligent transportation system ensuring optimal gap 
usage, e.g. based on intervehicle communication. Only in special cases, the task 
could be performed by a suitable traffic light control. 

In normal merging or intersection situations, there will always be coordination 
problems |16j when entering or crossing another flow, if the traffic volumes reach a 
certain level. This will cause inefficiencies in the usage of available road capacity, i.e. 
mixed flows will not be able to use the full capacity. Such effects can be modelled 
by specifying the corresponding permeabilities 7i(t) as a function of the merging 
flows, particularly the main flow or crossing flow. The deviation of "fi(t) from 1 
will then be a measure for the inefficiency. A particularly simple, phenomenological 
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specification would be 

^ = 1 + 0^-0.) ' < 62 > 
where the own outflow O2 supports a high permeability and the intersecting outflow 
0\ suppresses it. However, rather than using such a phcnomcnological approach, 
the permeability could also be calculated analytically, based on a model of gap 
statistics, since large enough vehicle gaps are needed to join or cross a flow. Such 
kinds of calculations have been carried out in Refs. [Pol 1501 I55"l 156] . 

4. Towards a Self-Organized Traffic Light Control. In Ref. [14], it has been 
pointed out that, for not too small arrival flows, an oscillatory service at intersections 
reaches higher intersection capacities and potentially shorter waiting times than a 
first-in- first-out service of arriving vehicles. This is due to the fact that the outflow 
of queued vehicles is more efficient than waiting for the arrival of other freely flowing 
vehicles, which have larger time gaps. For similar reasons, pedestrians are passing a 
bottleneck in an oscillatory way [2D] , and also two intersecting flows tend to organize 
themselves in an oscillatory way [TH [T7] . 

Therefore, using traffic lights at intersections is natural and useful, if operated 
in the right way. However, how to switch the traffic lights optimally? While this is 
a solvable problem for single traffic lights, the optimal coordination of many traffic 
lights [31] is a really hard (actually NP hard) problem [8]. Rather than solving a 
combinatorial optimization problem, here, we want to suggest a novel approach, 
which needs further elaboration in the future. The idea is to let the network flows 
self-organize themselves, based on suitable equations for the permeabilities -fi(t) as 
a function of the outflows Oi(t) and the numbers AiV,-(i) of delayed vehicles. 

Here, we will study the specification 

T^W = ,, n Mo a - 0l )-cD ( 63 ) 



T*(*)= i + oe^i-OaJ+cP ' (64) 



and 



which generalizes formula ([62]) . While the relative queue length 

D(t) = ANi (t) - AN 2 (t) (65) 

quantifies the pressure to increase the permeability 71 for road section 1, the outflow 
02(t) from the road section 2 resists this tendency, while the flow 0\(t) on road 
section 1 supports the permeability. The increasing pressure eventually changes the 
resistance threshold and the service priority. An analogous situation applies to the 
permeability 72 for road section 2, where the pressure corresponds to —D, which is 
again the difference in queue length, a, 6, and c are non-negative parameters, a may 
be set to 1, while c must be large enough to establish a sharp switching. Here, we 
have assumed c = 100. The parameter b allows to influence the switching frequency 
/, which is approximately proportional to b. We have adjusted the frequency / to 
the cycle time 

ycyc = 2t x ( 66 ) 

l-(A 1+ A 2 )/Q 

which results if the switching (setup) time ( "yellow traffic light" ) is t = 5s and a 
green light is terminated immediately after a queue has dissolved after lifting the 
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red light E The corresponding parameter value is 

500 



b = 



Q-(A 1 +A 2 ) 

Figure S] shows a simulation result for Ai/Q = 0.3 and A 2 /Q = 0.4. 



(67) 




100 150 



Figure 4. Illustration of the dynamics of self-organized oscilla- 
tions in the permeabilities and the resulting flows for a single inter- 
section with constant inflows (see text for details). Note that the 
road section with the higher inflow (arrival rate) is served longer, 
and its queues are shorter (see solid lines). 

The properties of the corresponding specification of the permeabilitities 7, (t) are 
as follows: 

• 7i(i) is non-negative and does not exceed the value 1. 

• For the sum of permeabilities and a > 1, we have 



7i + 72 = 



2 + a(e E + e~ E ) 



< 1 



l + a 2 + a(e E , . 
where we have introduced the abbreviation 

E = b(Oi - 2 ) + c(AJVi - AiV 2 ) . 
The sum is close to 1 for large absolute values of E, while for E 



(68) 



(69) 

the overall 
permeability 71 + 72 is small. 

For large enough values of ab and for c,A\,A 2 > 0, the equations for the 
permeability do not have a stable stationary solution. This can be concluded 
from 

dE _ h (d0 1 2 \ , fdAN! dAN 2 * 
~dt ~ ' 
together with 

= Ai-Oi(t) (71) 



dt dt 
dAN, 



V dt 



dt 



(70) 



dt 



3 If ATi and AT2 denote the green time periods for the intersecting flows 1 and 2, respectively, 
the corresponding red time periods for a periodic signal control are AT2 and ATi, to which the 
switching setup time of duration r must be added. From formula l|23| l and with Oi = Q we 
obtain ATi = (AT 2 + t)Q/(Q - Ai) and AT 2 = (ATi + t)Q/(Q - A 2 ). Using the definition 
T cyc _ ATi + t + AT 2 + r for the cycle time, we finally arrive at Eq. (1681 . 
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and 

Oi(t) = 7i (t) max[Q9(A^ > 0), A,} , (72) 

see Eqs. flU) and ©. As dD/dt = dANi/dt - dAN 2 /dt varies around zero, 
the same applies to D(t), which leads to oscillations of the permeabilities ji(t). 

• With the specification (|67[) of parameter b, the cycle time is approximately 
proportional to the overall inflow {A\ + ^2)- 

• The road section with the higher flow gets a longer green time period (see Fig. 

a. 

If the above self-organized traffic flows shall be transfered to a new principle of 
traffic light control, phases with 7j(t) ~ 1 could be interpreted as green phases and 
phases with 7i(£) ~ as red phases. Inefficient, intermediate switching time periods 
for certain choices of parameter values could be translated into periods of a yellow 
traffic light. 

5. Summary and Outlook. We have presented a simple model for conserved 
flows in networks. Although our specification has been illustrated for traffic flows 
in urban areas, similar models are useful for logistic and production system or even 
transport in biological cells or bodies. Our model considers propagation speeds of 
entities and congestion fronts, spill-over effects, and load-dependent transportation 
times. 

We have also formulated constraints for network nodes. These constraints con- 
tain several minimum and maximum functions, which implies a multitude of possible 
cases even for relatively simple intersections. It turns out that the arrival and de- 
parture flows of diverges have uniquely defined values, while merges or intersections 
have a set of feasible solutions. This means, the actual result may sensitively de- 
pend on the intersection design. For mathematical reasons, we have determined 
flow-optimizing solutions for two merging and two intersecting flows. However, it 
is questionable whether these solutions can be established in reality without the 
implementation of intelligent transport systems facilitating optimal gap usage: In 
many situations, coordination problems between vehicles in merging or intersection 
areas cause inefficiencies, which reduce their permeability. 

In fact, at not too small traffic volumes, it is better to have an oscillation between 
minimum and maximum permeability values. Therefore, we have been looking for a 
mechanism producing emergent oscillations between high and low values. According 
to our proposed specification (which is certainly only one of many possible ones) , 
the transition between high and low permeability was triggered, when the difference 
between the queue lengths of two traffic flows competing for the intersection capacity 
exceeded a certain value. The resulting oscillatory service could be used to define 
traffic phases. One potential advantage of such an approach would be that the 
corresponding traffic light control would be based on the self-organized dynamics 
of the system. Further work in this direction seems very promising. 
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